Experimental study of nonlinear dust acoustic solitary waves in a dusty plasma 
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The excitation and propagation of finite amplitude low frequency solitary waves are investigated 
in an Argon plasma impregnated with kaolin dust particles. A nonlinear longitudinal dust acoustic 
solitary wave is excited by pulse modulating the discharge voltage with a negative potential. It is 
found that the velocity of the solitary wave increases and the width decreases with the increase of 
the modulating voltage, but the product of the solitary wave amplitude and the square of the width 
remains nearly constant. The experimental findings are compared with analytic soliton solutions of 
a model Kortweg-de Vries equation. 

PACS numbers: 52.27.Lw, 52.27.Gr, 52.35. Sb 
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A dusty plasma, consisting of large macroparticles 
(usually the size of a few microns), immersed in an ion- 
ized gaseous medium with free floating electrons and 
ions, supports a rich variety of collective phenomena that 
have been the subject of a large number of investigations 
in recent times The massive dust particles, which 
also acquire a very high dynamic electric charge, bring 
about many significant changes in the overall collective 
behaviour of the system including the creation of new 
modes. The dust acoustic wave (DAW) is one such new 
mode [|J peculiar to this system and is an analog of the 
ion acoustic wave in a normal plasma with the dust par- 
ticles playing the role of the ions by providing inertia 
and the pressure contributions for sustaining the wave 
coming from both the lighter species (the electrons and 
the ions). The linear characteristics of DAWs have by 
now been well established both theoretically [3| and ex- 
perimentally [J]. There is also an enormous theoretical 
literature on the topic of nonlinear DAWs (e.g. on the 
formation of solitons, vortices, shock structures etc.)[^,0| 
but there has been relatively little experimental investi- 
gation 0, H, 0, [T3| of these nonlinear entities in a dusty 
plasma. Nakamura et al. Q showed that, an oscillatory 
ion-acoustic shock wave in a conventional two component 
Argon plasma transforms into a monotonic shock wave 
front when it propagates through a dusty plasma. Sam- 
sonov et al. [8| investigated the propagation and dissi- 
pation of longitudinal solitons in a 2D monolayer hexag- 
onal dust lattice of plastic monodisperse particles and 
invoked a linear chain theory, which included damping, 
dispersion and nonlinearity, to explain their experimental 
results. In a subsequent work the same group [9j reported 
the experimental observations of shock melting in a two 
dimensional dusty plasma and their experimental results 
were supported by molecular dynamics simulations. Pra- 
manik et al [l(| observed dust acoustic turbulence arising 
from nonlinear mode-couplings and associated harmon- 
ics generation. While there have theoretical predictions 
of the existence of DAW solitons in many space plasma 
contexts, to the best of our knowledge there has been no 
laboratory studies of DAW solitary structures in a reg- 
ular fluid state dusty plasma. In this Letter we present 
the first experimental evidence of such long lived struc- 



tures and describe their characteristics over a range of 
propagation velocities and spatial widths. We find that 
with the increase of the driving modulating voltage the 
solitary wave propagation velocity increases but its width 
decreases such that the product of the amplitude of the 
pulse and the square of its width remains nearly constant. 
We compare and interpret our results in terms of exist- 
ing theoretical soliton solutions of a model Korteweg- 
de Vries equation pertinent to our experimental parame- 
ter regime. 

The experiments have been carried out in a set-up 
shown in Fig. Q] which is similar to the set up described 
earlier in [TlJ. The stainless steel (SS) cylindrical cham- 
ber has eight radial ports and two axial ports for differ- 
ent purposes. The entire chamber is covered by a thin 
SS sheet to avoid stray arcing. Micron sized kaolin dust 
particles are sprinkled at the bottom of the SS cylindri- 
cal chamber. The vacuum chamber is initially pumped 
down to a base pressure of 10~ 3 mbar by a rotary pump. 
It is then purged with Argon gas using a precision nee- 
dle valve and subsequently pumped down again to base 
pressure. This process is repeated several times. Finally 
the pressure is kept constant at 1 mbar. A discharge 
is then produced between a rod shaped anode and the 
grounded vessel (used as a cathode) at an applied voltage 
of V a — 600 volts. The applied voltage is then reduced 
to 542 volts and the neutral gas pressure is gradually re- 
duced to 0.09 mbar to achieve the formation of a dense 
dust cloud. For such a discharge condition the average 
discharge current is usually ~ 105mA. 

The levitated dust particles are illuminated by a green 
Nd — Y AG diode laser. The laser light is spread into a 
sheet in the x — z plane by a cylindrical lens (kept hori- 
zontally) and forward scattered light from the dust cloud 
is used to visualize the dust particles. The scattered light 
from the dust particles are captured using a CCD cam- 
era (25 fps) which is kept at an angle of 15° to the y — z 
plane. Further, the distance measurement in the hori- 
zontal plane (x — z) is duly corrected to account for the 
geometrical effects arising from the small but finite angle 
of inclination of the CCD camera with respect to the z 
axis. Video frames are digitized by a frame grabber card 
with eight bit of intensity resolution and stored onto a 
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FIG. 1: A schematic of the experimental setup 



high speed computer. The average dust temperature (T^) 
is calculated from the random velocities of the particles 
by tracing single particle trajectories in different frames. 
The velocity, amplitude and width of the wave are all ob- 
tained from the still images. The plasma parameters like 
the ion density (m) and the electron temperature (T e ) are 
measured using a single Langmuir probe in the absence 
of dust particles inside a plasma. The electron density 
is estimated from the modified quasi-neutrality condition 
of the dusty plasma. A hot emissive probe is also used to 
measure the plasma potential and the floating potential. 
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FIG. 2: Typical voltage impulses used to excite dust acoustic 
solitary waves. 

Once the dust cloud is formed a short negative voltage 
pulse is superimposed over the discharge voltage to per- 
turb the equilibrium dust cloud. The dust cloud is then 
found to oscillate for the duration of the pulse length 
and to produce rarefaction and compressive structures 
corresponding to a dust acoustic wave. The oscillations 
quickly die down as soon as the pulse voltage stops. How- 
ever for a voltage pulse that exceeds a certain threshold 
value (60 volts in this case) the dust column at the end 
of its oscillation phase launches a small density perturba- 
tion. This newly born compressive structure (a dust den- 
sity perturbation) which constitutes a dust acoustic soli- 
tary wave (DASW) begins to move in the direction away 
from the anode and continues to move till the next pulse 
comes on (in this case about 1.5 to 2 sees later). The ap- 
plied pulse voltage width and the time between successive 
pulses for the excitation of DASWs is adjusted such that 



the solitary wave can propagate over a long distance. The 
typical values of the pulse width and inter-pulse period 
are 200 — 300 msec and 1.5 — 2.5 sec respectively and are 
shown in Fig. [5] 
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FIG. 3: The propagation of dust acoustic solitary waves, a) 
video image of a typical solitary structure and b) its pixel 
intensity profile, c) Snapshots of the time evolution of the 
solitary wave at intervals of 0.04 sec. 

Figure [3](a) , shows a typical solitary structure which 
propagates in the forward direction {z direction). Figure 
|3jb) shows the variations of its pixel intensity (averaged 
over x) along the z direction. This pixel intensity pro- 
file is obtained after carrying out a background reduction 
and is proportional to the particle density. Figure J3Jc), 
a typical time evolution plot of the solitary wave, shows 
that it moves with a constant speed for a long distance 
at a constant pulse height. These time sequence plots 
are made at intervals of 0.04 sec. From this figure, it is 
also clear that the amplitude (rid/ndo) and the full width 
at half maximum (FWHM) of a typical DASW remains 
nearly constant throughout its journey. The and ndo 
profiles shown in Fig. [3jc) correspond to pixel intensi- 
ties of DASWs and the equilibrium dust cloud respec- 
tively. The velocity of the DASWs is calculated from the 
distance traveled by a coherent structure between two 
consecutive frames. 

To further characterize the DASWs, we have varied 
the pulse height of the modulated discharge voltage by 
keeping the ON and OFF times of the pulse constant. 
The pulse height has been changed over a wide range 
from 60 to 120 volts. Figures HJa) andBJb) show typical 
DASW structures at two different peak voltages V pea k 
— 60 and 120 volts respectively. In Fig. Ufa), a com- 
pressive structure moves 5 mm in 0.28 sec. In the other 
figure (Fig.|4jb)), the structure moves 10 mm in the same 
time interval. It is clear from these figures that the dust 
acoustic solitary waves move with different velocities if 
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FIG. 4: Solitary structures at two different excited pulse volt- 
ages for (a) Vpp = 60 volts and (b) V pp = 120 volts. The wave 
propagates with higher velocity and amplitude but with lesser 
width when the applied voltage increases. 



the excitation voltages are different. 

We have also found that the amplitude of the com- 
pressive structure increases with the increase of the pulse 
height which is shown in Fig. IDJa) and Fig. 0Jb) . From 
these figures, it is also seen that the full width at half 
maximum (FWHM) decreases significantly with the in- 
crease of pulse height. 

In another set of experiments, we have observed the 
propagation of two solitary structures of different heights 
moving with different velocities. This happens when the 
driving pulse height exceeds a certain voltage (in this 
case it is 130 volts) and the resultant launched compres- 
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FIG. 5: The breaking of large amplitudes structures into two 
different small amplitudes which moves with different veloci- 
ties. 



sive pulse breaks up into two pulses. The amplitude of 
the leading pulse is seen to be a bit higher than the sec- 
ond one as shown in Fig. 02a). With time the distance 
between these two solitary waves is seen to keep on in- 
creasing indicating that the higher amplitude pulse is 
traveling faster (see figures Ufa) andEJb)). The above 
observations suggest that the dynamics of these solitary 
pulses are quite akin to theoretical results on compressive 
solitary waves of the Korteveg-de Vries kind. 

To carry out a comparison of our experimental results 
with such theoretical solutions we consider the follow- 
ing model K-dV equation that has been derived in the 
context of dust acoustic waves @. 
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where n d is the dust density perturbation normalized 
by the equilibrium dust density n<jo and the coefficients 
are, 



A = 
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Here £ and r are stretched coordinates such that £ = 
e x l 2 (x — Dot) and r = e 3 / 2 i with e representing a small- 
ness parameter measuring the weakness of the amplitude 
and vo is the wave phase speed normalized by the dust 
thermal velocity Co = (ZdhsTi/md) 1 ^ 2 ■ The laboratory 
frame time (t) and space (x) variables are normalized 
by the dust plasma frequency uj p d = (AimdoZ^e 2 /nid) 1 ^ 2 
and the plasma Debye length = (fcsTj/47rnoe 2 ) 1 / 2 . 
Zde = Qd and m,d are the dust charge and the mass of 
each dust grain respectively. Further S = nio/n e o and 
0i = Ti/Te where Ujo,Tj, (j = i,e) are the equilibrium 
densities and temperatures of the ions and electrons re- 
spectively. 

A solitary pulse (soliton) solution of Eqn.([T]) results 
from a balance between the nonlinear steepening induced 
by the convective term nP(drid/d^) and dispersive 

broadening due to the term d 3 n d / dt; 3 . Such a station- 
ary solution has the following exact analytic form, 
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where the amplitude n}^L — 38M/A and the width 
A s = y/4B/5M respectively. Here n — £ — 5Mt, with 
M = (vo/vda) = 1 + SM being the Mach number with 
respect to the linear dust acoustic velocity Vda- To es- 
timate the Mach number of the traveling pulse we have 
independently measured the linear phase velocity Vda of 
the DAW by modulating the DC discharge voltage with 
a very low amplitude sinusoidal voltage. Care was taken 
to ensure that the basic plasma conditions remained the 
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same as those for which we have studied the solitary wave 
propagation. The velocity of this linear dust acoustic 
wave is then calculated from consecutive frames in the 
interval of 0.04 sees and turns out to be ~ 1.67 cm/sec 
for our discharge conditions. It is to be noted that all 
the constants of the expression A (see eqn. ([2])) are posi- 
tive, making A to be positive. It indicates that the dusty 
plasma supports dust acoustic solitary waves with a neg- 
ative potential as 4>m = ~ v o n< dm- 
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FIG. 6: Variation of a) solitary amplitude, b) solitary width 
and c) product of solitary amplitude and the square of its 
width with Mach number. Solid line represents the solution of 
K-dV equation and open circles are the experimental points. 

The width and amplitude of a K-dV soliton solution 
([3]) have a definite dependence on the pulse propagation 
speed. These theoretical dependencies (solid line) of am- 
plitude, width and the product of the square of the width 



and the amplitude with the Mach number are plotted in 
Fig. ©(a)-(c) along with experimental results (empty 
circles). All the constants were estimated using experi- 
mentally measured parameters like Z d = 3 x 10 3 , T e = 8 
eV, n % = 7x 10 13 m" 3 , n d = 1 x 10 10 itT 3 , m d = 1 x lO" 13 
kg. The ion temperature in our experiments is bit high 
(=0.3 eV) but consistent with similar earlier measure- 
ments reported by Thompson et al. As can be 
seen, the experimental results agree quite well with the 
theoretical estimates particularly for low amplitudes and 
low Mach numbers. The deviations at higher amplitudes 
(and Mach numbers) can be attributed to several factors. 
First of all the K-dV model itself ceases to be valid be- 
yond a certain amplitude and corrections due to higher 
order nonlinearities become important. Among other 
physical factors that could be contributing to the devia- 
tion one should comsider dust charge fluctuation effects 
as well as strong coupling effects, as have been pointed 
out in some past theoretical studies [l3[. However a re- 
liable estimate of these effects is difficult to quantify at 
this time because of experimental limitations (in making 
precise measurements of the involved quantities) as well 
as due to lack of reliable theoretical estimates of these 
effects in the nonlinear regime. 

To conclude, we report experimental observations of 
long lived low frequency solitary wave pulses in a labo- 
ratory dusty plasma. These pulses that correspond to 
nonlinear dust acoustic waves, are externally excited in a 
controlled manner and are investigated over a wide range 
of amplitudes and propagation velocities. Their dynam- 
ical properties are found to agree quite well, particularly 
at low amplitudes and low Mach numbers, with those of 
solitonic solutions of a model K-dV equation describing 
the propagation of nonlinear dust acoustic waves. The 
observed deviations at larger amplitudes indicate the in- 
fluence of higher nonlinearities and/or the presence of 
additional physical effects that merit future experimen- 
tal investigations. 



[*] Corresponding author: E-mail address: banerjee_pintu20 
02@yahoo.com (P. Bandyopadhyay) 

[1] F. Verheest, Waves in Dusty Space Plasmas, Kluwer Aca- 
demic, Dordrecht, 2000; V.E. Fortov, A.G. Khrapak, 
S.A. Khrapak, V.I. Molotkov, O.F. Petrov, Phys. Usp. 
47 (2004) 447; and references therein. 

[2] N. N. Rao, P. K. Shukla and M. Y. Yu, Planet. Space 
ScL, 38, 543 (1990). 

[3] P.K. Shukla, A. A. Mamun, Introduction to Dusty Plasma 
Physics, Institute of Physics, Bristol, 2002; and refer- 
ences therein. 

[4] R. L. Merlino, A. Barkan, C. Thompson and N. 
D'Angelo, Plasma Phys. Control. Fusion, 39, A421 
(1997); 

[5] N. N. Rao, Physica Scripta, T75, 1179 (1998). 

[6] A. V. Volosevich and C. V. Meister, Contr. to Plasma 
Phys., 42, 61 (2002); P. K. Shukla and A. A. Mamun, 
New J. Phys., 5, 17 (2003); P. K. Shukla, Phys. Plas- 



mas, 10 1619 (2003); S. I. Popel, A. G. Golub and T. V. 

Losseva et al, Phys. Rev. E, 67, 056402 (2003) 
[7] Y. Nakamura, H. Bailung and P. K. Shukla, Phys. Rev. 

Lett. 83, 1602 (1999). 
[8] D. Samsonov, A. V. Ivlev, R. A. Quinn and G. Morrill, 

Phys. Rev. Lett. 88, 095004 (2002). 
[9] D. Samsonov, S. K. Zhdanov, R.A. Quinn, S. I. Popel 

and G. E. Morfill, Phys. Rev. Lett., 92, 255004 (2004). 
[10] J. Pramanik, B. M. Veeresha, G. Prasad, A. Sen and P. 

K. Kaw, Phys. Lett. A, 312, 84 (2003). 
[11] P. Bandyopadhyay, G. Prasad, A. Sen and P. K. Kaw, 

Phys. Lett. A, 368, 491 (2007). 
[12] C Thompson, A Barkan, N D'Angelo and R. L. Merlino, 

Phys. Plasmas, 4, 2331 (1997). 
[13] S. Ghosh et al, Phys. Plasmas 9, 1150 (2002); P.K. Kaw 

and A. Sen, Phys. Plasmas 5, 3552 (1998). 



